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Abstract 

We calculate the contribution of the higher twist Feynman diagrams to the large-py inclusive pion 
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higher twist subprocess cross section and it's resummed expression are found. We compare the 
resummed higher twist cross sections with the ones obtained in the framework of the frozen coupling 
approximation and leading twist cross section. We discuss the phenomenological consequences 
of possible higher-twist contributions to the pion production in proton-proton collisions within 
holographic QCD. 
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I. INTRODUCTION 



One of the most interesting new developments in hadron physics is the application of the 
AdS/CFT correspondence l| to nonperturbative QCD problems 2|-|8|. It is well know that 
AdS/CFT gives an important insight into the viscosity and other global properties of the 
hadronic system formed in heavy ion collisions [9]. The essential ansatz for the application of 
AdS/CFT to hadron physics is the indication that the QCD coupling a s (Q 2 ) becomes large. 
Therefore conformal symmetry can be applied to, for example solutions of the QCD Dyson 
Schwinger equations and phenomenological studies of QCD couplings based on physical 
observables such as r decay and the Bjorken sum rule show that the QCD function 
vanishes and a s (Q 2 ) become constant at small virtuality; i.e., effective charges develop an 
infrared fixed point. Fully exploiting the gauge/gravity correspondence to produce a model 
for real strong interaction physics- a method called " holographic QCD" or " AdS/QCD"-may 
be attempted either through a top-dawn approach starting with a particular string theory 
and choosing a background that naturally produces QCD-like properties, or a bottom-up 
approach starting with real QCD properties and using them to obtain constraints on viable 
dual gravity theories. 

The first attempts have been made for constructing phenomenological holographic models 
of QCD [4]. Surprisingly simple models consisting of gauge theory in an anti-de-Sitter space 
interval have turned out to provide a remarkably good description of the meson sector of 
QCD. Therefore, it will be interesting that the calculation of the higher twist effects within 
holographic QCD in proton-proton collisions in the running coupling approach. 

The large-order behavior of a perturbative expansion in gauge theories is inevitably dom- 



inated by the factorial growth of renormalon diagrams |10Ml3|. In the case of quantum 
chromodynamics (QCD), the coefficients of perturbative expansions in the QCD coupling 
a s can increase dramatically even at low orders. This fact, together with the apparent free- 
dom in the choice of renormalization scheme and renormalization scales, limits the predictive 
power of perturbative calculations, even in applications involving large momentum transfer, 
where a s is effectively small. 

In this work we apply the running coupling approach 14| in order to compute the effects 
of the infrared renormalons on the pion production in proton-proton collisions within holo- 



graphic QCD. This approach was also employed previously 15Ml8| to calculate the inclusive 
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meson production in proton-proton and photon-photon collisions. 

For the calculations of the higher-twist cross sections on the dependence of wave func- 
tions of pion, we used the holographic QCD prediction & ho i and pion asymptotic 
wave functions Q asy (x) [22|] from the perturbative QCD evolution. Theoretically obtained 
predicted results within holographic QCD were compared with results of the perturbative 
QCD which obtained by the running coupling and frozen coupling constants approaches. 

The frozen coupling constant approach in Refs. was used for calculation of inte- 

grals, such as 

1 — x 

It should be noted that, in pQCD calculations, the argument of the QCD coupling constant 
(or the renormalization and factorization scale) Q 2 should be taken equal to the square of 
the momentum transfer of a hard gluon in a corresponding Feynman diagram. But definition 
of a s (Q 2 ) suffers from infrared singularities. Therefore, in the soft regions as X\ — > 0, and 
£2 —> 0, integrals (1.1) diverge and for their calculation some regularization methods are 
needed for a s (Q 2 ) in these regions. Investigation of the infrared renormalon effects in various 
inclusive and exclusive processes is one of the most important and interesting problem in 
the perturbative QCD. It is known that infrared renormalons are responsible for factorial 
growth of coefficients in perturbative series for the physical quantities. But, these divergent 
series can be resummed by means of the Borel transformation [lo| and the principal value 
prescription {27]. Studies of higher- twist and renormalon effects also opened new prospects 
for evaluation of power suppressed corrections to processes characteristics. 

We organize the paper as the follows. In Section [Til we provide some formulas for the 
calculation of the contributions of the higher twist and leading twist diagrams. In Section 
Hilt we present the formulas and analysis of the higher-twist effects on the dependence of 
the pion wave function by the running coupling constant approach, and in Section IIVI the 
numerical results for the cross section and discuss the dependence of the cross section on the 
pion wave functions are presented. Finally, some concluding remark are stated in Section PVT. 



3 



II. HIGHER TWIST AND LEADING TWIST CONTRIBUTIONS TO INCLU- 
SIVE REACTIONS 



The higher- twist Feynman diagrams, which describe the subprocess q\ + q 2 — > 7i + (tt~) +7 
for the pion production in the proton-proton collision are shown in Fig.l. Th e amplitude 
for this subprocess can be found by means of the Brodsky-Lepage formula 



28J 



M(s, i) 



(2-1) 



dx x \ dx 2 5(l —Xi - x 2 )$ n (x lj x 2 ,Q 2 )T H (s,i;x ll x 2 ). 
'0 Jo 

In Eq.(2.1), T H is the sum of the graphs contributing to the hard-scattering part of the 
subprocess. 

The Mandelstam invariant variables for subprocesses q\ + q 2 — > 7i + (tt~) + 7 are defined 



as 



s = ( Pl + P2 ) 2 , i=(p 1 -p 7I ) 2 , u = {p 1 -p~ t ) 2 . 



(2.2) 



The pion wave functions predicted by AdS/QCD [19h21| and the PQCD evolution [22j has 
the form: 



*JS(X) = —f K y/ X (l- X ), 



y/x(l -x)exp y- 2 kjx{l-x) J ' = ^-M 1 ~ x ) ( 2 - 3 ) 



®vsbgl{ x ) 

where f n is the pion decay constant. 

The cross section for the higher-twist subprocess q\q 2 —> 7r + (71-^)7 is given by the expres- 
sion 

^,i, fi) = !5!^M![l + Il (2 .4) 
dV ! 27 s 3 [v? t 2 \ V ' 

where 



D(i, u) = e\t I dx 



a a {Ql)^{x ,Q\ 
1 — x 



2\1 



1 

+ e 2 u I dx 
'0 



g s (Ql)%(x,Q 

1 — X 



2\1 



(2.5) 



In the Eq.(2.5) Q\ = (x — l)u and Q\ = —xi represent the momentum squared carried 



by the hard gluon in Fig.l, t\{e 2 ) is the charge of qi(q 2 ) and Cp 



The higher-twist 



contri 
is 



Dution to the large-pr pion production cross section in the process pp — > -k + (-k )+7+X 
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s m T = E 



da 
d 3 p 



dxidx 2 G qi / hl (xi)G m / h2 (x 2 )--^(qq -> ir>y)5(s + i + u). (2.6) 



J 



7T dt 
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We denote the higher-twist cross section obtained using the frozen coupling constant ap- 
proach by (£^ T )°. 

Regarding the higher-twist corrections to the pion production cross section, a compar- 
ison of our results with leading-twist contributions is crucial. We take two leading-twist 
subprocesses for the pion production:(l) quark-antiquark annihilation qq — > #7, in which 
the g — > 7r + (7r~) and (2) quark-gluon fusion, qq — > qj, with subsequent fragmentation of the 



final quark into a meson, q — > 7r + (7r ) 15|, [17]. 



III. THE HIGHER TWIST MECHANISM IN HOLOGRAPHIC QCD AND IN- 
FRARED RENORMALONS 

The main problem in our investigation is the calculation of integral in (2.5) by the running 
coupling constant approach within holographic QCD and also discussion of the problem of 
normalization of the higher twist process cross section in the context of the same approach. 
Therefore, it is worth noting that, the renormalization scale (argument of a s ) according to 
Fig.l should be chosen equal to Q\ = (x — l)u, Q\ = —xt. The integral in Eq.(2.5) in the 
framework of the running coupling approach takes the form 

j, 2 v = f 1 a s( A /4 )$M(z,/4)^ ^ 3 - 

R ° Jo 1 — X 

The a s (\fi 2 Ro ) has the infrared singularity at x — > 1, for A = 1 — x or x — > 0, for A = x and 
so the integral (3.1) diverges. For the regularization of the integral, we express the running 
coupling at scaling variable a s (Xfi 2 Ro ) with the aid of the renormalization group equation 



in terms of the fixed one a s (Q 2 ). The so 
running coupling a = a s /ir has the form 

a(A) 



ution of renormalization group equation for the 
27] 



a 



1 + a— In A 
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(3.2) 



Then, for a s (\Q 2 ), we get 



where t = 4n/a s (Q 2 )(3 = 4/a/3 . 

Having inserted Eq.(3.3) into Eq.(2.5) we obtain 

D{i , & ) = e lt " f dx ".^)^,Q\) + e2fi f /.(AA)*M(x.Cj) 
Jo 1 - x Jo l-x 

5 



eitasi-ujtx / dx- 

Jo (1 -x){ti + lnA) 



+ e 2 Ma s (- t)t 2 / g^t— 



(3.4) 



-x)(* 2 + lnA) 
where £i = Att /a s {— u)[3 and i 2 = 47r/a< j (— f)/3 . 

Although the integral (3.4) is still divergent, it is recast into a suitable form for calculation. 
Making the change of variable as z = In A, we obtain 



D(t,u) = eita s (—u)ti / dx 



+ e 2 wa s (— f)t 2 / dx 



$m(x,QI 



/ (l-x)(h + z) ' sy ~'"J Q (l-x)(t 2 + z) 
In order to calculate (3.5) we will apply the integral representation of l/(t + z 



(3.5) 
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32]. 



1 



(t + z) 



- {t+z)u du, 



(3.6) 



gives 



Z 1 f°° ®Jx,Ql)e-^ +z >dudx . . . f l f°° <$>Jx,Q 2 2 )e-^ +z >dudx 
D(t,u) = eita s (-«)ti / / v '^y r +e 2 wa s (-i)t 2 ' ' 



JO 



1-x 



JO 



1-x 
(3.7) 



In the case $^(x) for the D(i,u) it is written as 



L>(£, w) 



2' 2 



3 1 



and for ^sy( x ) wave function 



Po JO 



1 — it 2 — it 



4\/37r/ 7r 



e 2 -u 



/io 



-t 2 « 



1 



1 



1 — u 2 — it 
(3.9) 



where B(a,/3) is Beta function. The structure of the infrared renormalon poles in Eq.(3.8) 
and Eq.(3.9) strongly depend on the wave functions of the pion. To remove them from 
Eq.(3.8) and Eq.(3.9) we adopt the principal value prescription. We denote the higher-twist 
cross section obtained using the running coupling constant approach by (S^ T ) res . 

IV. NUMERICAL RESULTS AND DISCUSSION 



In this section, we discuss the higher-twist contributions calculated in the context of 
the running and frozen coupling constant approaches on the dependence of the chosen pion 
wave functions in the process pp — > 7r + (or 7r~)7 + X. In numerical calculations for the quark 
distribution function inside the proton, the MSTW distribution function [331 ] . 



and quark fragmentation 



and the gluon 



341 ] functions into a pion have been used. The results of our 



numerical calculations are displayed in Figs. 2-14. Firstly, it is very interesting comparing 
the higher-twist cross sections obtained within holographic QCD with the ones obtained 
within perturbative QCD. In Fig. 2 and Fig. 3 we show the dependence of higher-twist cross 
sections (E^ 1 ) , (Y,^+) res calculated in the context of the frozen and running coupling 
constant approaches as a function of the pion transverse momentum p? for different pion 
wave functions at y — 0. It is seen from Fig. 2 and Fig. 3 that the higher-twist cross section is 
monotonically decreasing with an increase in the transverse momentum of the pion. In Fig. 4- 
Fig.7, we show the dependence of the ratios / \Y? HT ) , (Ef T ) re 7(Ef+ T )°, (Ef + T )°/(E^) 
and (E^ r ) res / as a function of the pion transverse momentum p T for $^ oZ (x), 3>£(x) and 
&vsbgl( x ) pi° n wave functions. Here E^+ is the leading-twist cross section, respectively. As 
shown in Fig. 4, in the region 2 GeV/c < pt < 30 GeV/c resummed higher-twist cross section 
for (&^°\x) is suppress by about half orders of magnitude relative to the resummed higher- 
twist cross section for also higher- twist cross section for 3>£(x) is suppress by about 
half orders of magnitude relative to the higher-twist cross section for &v°sbgl( x ) ■ 111 Fig- 5 
and Fig.6, we show the dependence of the ratios (Ef T ) re 7(Ef+ T )°, and (Ef+ T )°/(E^) as a 
function of the meson transverse momentum pr for the ^ ol (x), <&%(x) and &v°sbgl( x ) pion 
wave functions. It is observed from Fig. 5 that, the ratios for all wave functions decrase with 
an increase in the transverse momentum of pion, but in the region 2 GeV/ c < pr < 8 GeV/ c 
resummed higher-twist cross section for is suppress by about 2-3 and 2-4 orders of 
magnitude relative to the higher-twist cross sections calculated in the context of the frozen 
coupling method for and &v°sbgl( x )i respectively. Also, resummed higher- twist cross 
section for &%(x) is suppress by about 1-3 orders of magnitude relative to the higher- twist 
cross sections calculated in the context of the frozen coupling method for *&vsbgl( x )- But, 
as shown in Fig.6, the ratios decrease with increasing in the pt transverse momentum of the 
pion and has a minimum approximately at the point pt = lSGeV/ c, then the ratios increase 
with increasing in the pt transverse momentum of the pion. In Fig. 7 is similar to Fig.6 
with an exception that there in a minimum approximately at the point pr = 25GeV/c. In 
Fig.8 - Fig.10, we have depicted higher-twist cross sections (E^ 1 ), and ratio (S^ ) y)/(E^-/ r ), 
as a function of the rapidity y of the pion at y/s = 62.4 GeV and pt = 4.9 GeV jc. Figures 
show that higher-twist cross section and ratios have a different distinctive. The resummed 
higher-twist cross section for Q^ ol (x) has a maximum approximately at the point y = —1.92. 
However in this point higher-twist cross section for ®vsbgl( x ) nas a minimum. As shown in 
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Fig.9 in the region (-2.52 < y < -1.92) the ratios {Y^ l T f / {Y? HT f and (^vsbgl)° /(^ht)° 
increase with an increase of the y rapidity of the pion and has a maximum approximately 
at the point y = 1.22, but ratio (S^j,) '/ '(^v°sbgl)° nas a minimum in this point. As is seen 
from Fig. 10 ratios R = (S^ r ) res /(S^ r )°, for all wave functions increase with an increase 
of the y rapidity of the pion and has a maximum approximately at the point y = —1.92. In 
the region —2.52 < y < —1.92 resummed higher-twist cross section is suppress by about half 
order of magnitude relative to the higher-twist cross section calculated in the framework of 
the frozen coupling approach. Besides that, the ratio decreases with an increase in the y 
rapidity of the pion. As is seen from Fig. 10, the ratio R depends on the choice of the pion 
wave function. Analysis of our calculations shows that (E^) , (S^ r ) res higher-twist cross 
sections and ratio sensitive to pion wave functions predicted holographic and perturbative 
QCD. 

We have also carried out comparative calculations in the center-of-mass energy y/s = 
200 GeV and obtained results are displayed in Fig.ll-Fig.14. Analysis of our calculations at 
the center-of-mass energies y/s = 62.4 GeV and y/s = 200 GeV, show that with increasing in 
the beam energy contributions of higher-twist effects to the cross section decrease by about 
one-two order. As is seen from Fig. 5, Fig.7, Fig.10 and Fig.13 that infrared renormalon 
effects enhance the perturbative predictions for the pion production cross section in the 
proton-proton collisions about 2-3 order. This feature of infrared renormalons may help the 
explain theoretical interpretations with future experimental data for the pion production 
cross section in the proton-proton collisions. In our calculations, the higher-twist cross 
section of the process the dependence of the transverse momentum of pion appears in the 
range of (10~ 8 -j- 10~ 26 )mfe/ GeV 2 . Therefore, higher-twist cross section obtained in our work 
should be observable at RHIC. 

V. CONCLUSIONS 

Proton-proton collisions are known are known to be the most elementary interactions 
and form the very basis of our knowledge about the nature of high energy collisions in 
general. Physicists, by and large, hold the view quite firmly that the perturbative quantum- 
chromodynamics provides a general framework for the studies on high energy particle- 
particle collisions. Obviously, the unprecedented high energies attained at Large Hadron 
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Collider offer new window and opportunities to test the proposed QCD dynamics with its 
pros and cons. However, we should remember that LHC opens a new kinematical regime 
ah high energy, where several questions related to the description of the high-energy regime 
of the QCD. Consequently, studies of proton-proton interactions at the RHIC and LHC 
could provide valuable information on the QCD dynamics at high energies. In this work the 
single meson inclusive production via higher twist mechanism within holographic QCD are 
calculated . For calculation of the cross section, the running coupling constant approach is 
applied and infrared renormalon poles in the cross section expression are revealed. Infrared 
renormalon induced divergences is regularized by means of the principal value prescripton 
and the resummed expression (the Borel sum) for the higher-twist cross section is found. 
It is observed that, the resummed higher-twist cross section differs considerably from that 
found using the frozen coupling approximation, in some region. The following results can be 
concluded from the experiments; the higher-twist contributions to single meson production 
cross section in the proton-proton collisions have important phenomenological consequences, 
the higher-twist pion production cross section in the proton-proton collisions depends on the 
form of the pion model wave functions and may be used for their study. Also the contribu- 
tions of renormalons effects within holograpich QCD in this process are essential and may 
help to analyse experimental results. Further investigations are needed in order to clarify 
the role of higher-twist effects in QCD. Especially, the forthcoming RHIC and LHC mea- 
surements will provide further tests of the dynamics of large-pr hadron production beyond 
the leading twist. 
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FIG. 1: Feynman diagrams for the higher twist subprocess, qiq2 — > vr + (or ir ) r y. 
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FIG. 2: Higher- twist 7r + production cross section (£-^p )° as a function of the pr transverse mo- 
mentum of the pion at the cm. energy \fs = 62.4 GeV . 
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FIG. 3: Higher-twist 7r + production cross section (Yl H +) res as a function of the pt transverse 
momentum of the pion at the cm. energy a/s = 62.4 GeV. 
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FIG. 4: Ratio (Y,^ T ) / (Y? HT ) , where higher-twist contribution are calculated for the pion rapidity 
y = at the cm. energy y/s = 62.4 GeV as a function of the pion transverse momentum, px- 
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FIG. 5: Ratio (E +) res / (E +)° , as a function of the px transverse momentum of the pion at the 
cm. energy y/s = 62.4 GeV. 
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FIG. 6: Ratio (S^) /(£+), where higher-twist contribution are calculated for the pion rapidity 
y = at the cm. energy ^fs = 62.4 GeV as a function of the pion transverse momentum, px- 
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FIG. 7: Ratio (S^ r ) res /(S^+ )°, as a function of the px transverse momentum of the pion at the 
cm. energy ^fs = 62.4 GeV . 
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FIG. 8: Higher- twist 7r + production cross section (T,^+) , as a function of the y rapidity of the 
pion at the transverse momentum of the pion pt = 4.9 GeV/c, at the cm. energy y^s = 62.4 GeV. 




FIG. 9: Ratio (S^) /(S^ r ) , as a function of the y rapidity of the pion at the transverse 
momentum of the pion pt = 4.9 GeV/c, at the cm. energy y/s = 62.4 GeV. 
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FIG. 10: Ratio (£^+ , ) res /(^viO°> as a function of the y rapidity of the pion at the transverse 
momentum of the pion p? = 4.9 GeV/c, at the cm. energy \fs = 62.4 GeV. 
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FIG. 11: Higher-twist ir + production cross section {Jj H ^) res as a function of the pt transverse 
momentum of the pion at the cm. energy y/s = 200 GeV. 
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FIG. 12: Ratio (S^)°/(E^ r )°) as a function of the px transverse momentum of the pion at the 
c.m.energy \fs = 200 GeV. 
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FIG. 13: Ratio (S^) res /(S^)°, as a function of the pt transverse momentum of the pion at the 
c.m.energy yfs = 200 GeV . 
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FIG. 14: Ratio (T,^+) res /(S^?") , as a function of the y rapidity of the pion at the transverse 
momentum of the pion pt = 15.5 GeV/c, at the cm. energy y/s = 200 GeV. 
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